A texture tensor to quantify deformations 
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Under mechanical deformation, most materials exhibit 
both elastic and fluid (or plastic) responses. No existing for- 
malism derived from microscopic principles encompasses both 
their fluid-like and solid-like aspects. We define the statisti- 
cal texture tensor to quantify the intuitive notion of stored 
deformation. This tensor links microscopic and macroscopic 
descriptions of the material, and extends the definition of elas- 
tic strain. 

A typical mechanical experiment applies a given 
macroscopic distortion to a test sample and measures 
the resulting macroscopic force exerted by the material, 
or vice versa. The goal is to find the constitutive equa- 
tion, which relates the macroscopic stress tensor to an 
independent descriptor of the material's response. We 
know the relevant descriptor to use for two extreme cases, 
elastic solids and isotropic fluids: the gradient of the dis- 
placement and the velocity field, respectively. 

While stress is unambiguously defined [1,2], strain ad- 
mits more than one definition. Classical linear (or even 
non-linear) elasticity operationally defines the strain by 
comparing the current microscopic state to a fixed micro- 
scopic reference state [3-5]. But most materials, having 
both an internal structure which stores elastic energy, 
and the flexibility to allow rearrangements, lie between 
ideal fluids and purely elastic solids. Not only do we lack 
their exact constitutive relations, we do not know what 
descriptors apply. 

In this paper we propose an operational definition of 
the deformation which we can measure in experiments 
and simulations in terms of averages of microscopic quan- 
tities: the statistical texture tensor, a state function of the 
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material. Its variations measure the elastic strain of an 
object under arbitrary deformations, without requiring 
the microscopic details of a reference state. 

Contrary to Ref. [6] which starts from a coarse grained 
mass density, which is appropriate for granular materials, 
we consider a network structure. Our generic material is 
a network of sites connected by links, which can detach 
from and reattach to other sites (Fig. 1). Site and link 
definitions depend on the material: 

(i) In cellular patterns, a site is the meeting point 
of cells. Two sites connect if their cells share an edge. 
This concerns for instance grain boundaries in crystals 
[7] , compact 2D or 3D aggregates of biological cells, or 
Voronoi tesselation. 

(ii) If links are physical objects, sites need not be. In 
liquid foams, links are bubble edges, while sites are the 
vertices where they meet. In gels of polymers, links are 
macromolecules; sites are knots. 

(iii) If sites are undeformable objects, as in hard gran- 
ular materials, they link if their separation is less than a 
cut-off distance. Since such cut-off is arbitrary, it must 
be chosen consistently throughout the analysis. 

(iv) When two sites exert a force on each other, they 
are linked. This is the case e.g. for atoms or molecules in 
crystalline and amorphous solids, or deformable granular 
materials. We then have to specify an cut-off on their 
interaction force (Cauchy scheme): again, such arbitrary 
cut-off must be chosen consistently. 

Note that, even in cases (ii) and (iv), the strain re- 
mains purely geometric and does not depend explicitly 
on stresses and forces. This is seen in 3D foams where 
the films determine the stress [8] while the edges deter- 
mine the statistical strain. 

For simplicity, we consider only point-like isotropic 
sites, i.e. the average link length is much larger than 
the site size, and the unstressed material is mesoscopi- 
cally isotropic [9]. Both restrictions are inessential: we 
could relax them by extending our definitions. 

At the microscopic level, we describe our material by 
the positions {r s } of the sites {s} in a d-dimensional 
space (Fig. 1). The topology T is the list of all pairs (s, s') 
of sites connected by a link I = (£ x , £ y , . . . ) = ry — r s . The 
topology changes when perturbations create or destroy a 
link or a site. 
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FIG. 1. A network of interconnected sites. The represen- 
tative volume element V at position R is a square in two di- 
mensions, its volume is V. Microscopically, r is the position 
of a site, and there are here three links (unoriented vectors I) 
per site. 

A continuous description interpolates between this mi- 
croscopic level of description and a macroscopic descrip- 
tion in terms of the shape and size of the network bound- 
aries. We define an intermediate mesoscopic level (Fig. 
1). We cut the network into representative volume ele- 
ments, V(R) of volume V, with positions R: small enough 
that their properties are constant over the box; but large 
enough that each piece contains enough links to compute 
statistical properties and average out microscopic details. 
In each volume element we define all statistical quanti- 
ties as averages over all links in V(R). The number of 
neighbours of each site must be much smaller than the 
number of links in V(R). We thus require that the topol- 
ogy remains "short-range:" while the neighbours s' of a 
site s may change, their number must remain bounded. 

Our fundamental definition is the texture tensor : 



M (R) 



(i) 



Here (£<E) T)^ = i^j is the "dyadic" (or "outer" or "ten- 
sor") product. M is symmetric, and has positive eigen- 
values. M occurs in many different physical contexts: 
the Sterner Tensor [10-12], order parameters of nemat- 
ics [13], molecular moments of inertia ( [14] p. 116-119), 
or textures of granular materials [15]. The texture tensor 
quantifies our mental image of a network; large eigenval- 
ues correspond to directions of stretching (Fig. 2). 

This definition (1) encompasses all materials where a 
mesoscopic scale exists, i.e., most cases. It requires that a 
thermodynamic limit exists for all extensive and intensive 
quantities, i.e. statistics on larger volumes decrease the 
relative amplitude of fluctuations (these microscopic fluc- 
tuations might remain visible on the scale of the sample 



size if the sample is small [16] or if the lattice is ordered 
[17]). We assume no correlations between the volume 
elements, which we must check. Materials that display 
networks of forces, avalanches or fractures thus require 
careful treatment. For materials in crgodic steady flow, 
time averages also reduce fluctuations and the network 
acts as a continuous medium down to scales as small as 
the average link length (Fig. 2). 




FIG. 2. Texture tensor M (eq. 1) superimposed on a snap- 
shot of a simulated two-dimensional foam. The foam flows 
steadily from left to right around a fixed round obstacle. We 
calculate the texture tensor over a box almost as small as the 
bubble size, by averaging over 50 successive images. The tex- 
ture tensor is symmetric with two positive eigenvalues. We 
represent it as an ellipse, with the long axis (resp. small 
axis) proportional in length to the largest (resp. smallest) 
eigenvalue and pointing along the corresponding eigenvector. 
We see compression defomation in front of the obstacle, and 
stretching behind it. 

If each region has an isotropic reference state [9] , then 

Mq= Mq Id, where / is the identity tensor in d di- 
mensions. All its relevant information lies in the scalar 
Mo = (^o) /d, i-e. in the mean squared link length in 
this state [18]. 

We now define the statistical strain U as: 



= {R) s log M — log Mq ^ 



(2) 



where the tensor log M has the same axes as M but 
eigenvalues equal to the logarithm of those of M. 

Eq. (2) is operational: U quantifies the deformations 
visible locally in each region of an image such as exten- 
sion, compression, shear, or dilation (Fig. 2). U also sat- 
isfies other requirements for the definition of strain. It 
is invariant under rotation, translation and transposition 
of indices. It is a true function of state: it depends ex- 
plicitly only on the current state and the reference state, 
not on the detailed history of the material. It is meso- 
scopic: it reflects only statistically significant features of 
the network. Its definition (eq. 2) accomodates topolog- 
ical changes and thus the plastic regime. 
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Finally, we prove that whenever the classical defini- 
tion of the strain tensor applies, it coincides with our 

statistical strain U. Under an arbitrary elastic deforma- 
tion, from an initial, equilibrium, isotropic state to a final 
state, we have: 

U= X - (fi + w) + O (w 2 ) , (3) 

Here we do not make any assumption about the micro- 
scopic displacements. The displacement field u(R) is the 
difference between states of the region V(R) initially at 
position R. When u(R) is differentiable, the distortion 
tensor w= [V<S>«]* = dui/drj quantifies the deformation 
from initial to final state. The proof goes as follows. 

As a first step, we have to prove that the average over 
the surface element n dS (or line elements, if in two di- 
mension), denoted (.}", of any function g(£) of the links 
which actually cross this surface, relates to its bulk aver- 
age, denoted (.), as: 

Vn, (gf (£ j )n j = (g£ j )n j . (4) 

Here (£j } depends on the orientation chosen for the link: 
it must be chosen consistently along the surface, for in- 
stance along the orientation of n. On the opposite, if 
g is an odd function of the link I, i.e. g{—£) = —g{£), 
then the quantity g{£)£ does not depend on the arbitrary 
orientation choosen for the link; in that case, (g£j) does 

not depend on n, so that we can omit n in what follows. 
To prove this first step, eq. (4), we remark that for 

any elementary section dS=n dS of C, oriented by the 

unit vector n, the probability for a link to cross dS is 
proportional to: 

pi n dSP(£)d 3 £, 

where p is the average density of links, and P is the 
probability distribution function of the links £. Hence the 

sum of g taken over all links crossing dS is: g(£) = 

dS 

p fff g{£) IndS P{£) d 3 £. In particular, with g = 1 we 
obtain the number of links crossing dS: 




I. n dSP(£)d 3 £. 



Combining both proves eq. (4). 

At this point, it is interesting to apply eq. (4) to the 
forces. By definition [3,4], the stress a is the tensor such 
that, V n: 

(TijUj ds — y~v t , 

dS 



where r is the tension (force) supported by the link £. 
With g{£) = n, eq. (4) yields: 

~a= p{7 ®/), 

providing a statistical interpretation of the stress tensor 
as an average over all links [8] . 

Now, we can achieve the proof of eq. (3). By definition 
of u and w, we can write: 

(6^ = (5r) n (R +(/)") - (<5r)" (r) 
= u(R+(^)-u(R) 
= w it)" +0(w 2 ). (5) 

With g{£) = Sii, eq. (4) yields: 

which inserted in eqs. (1) and (5) yields: 

M = M + M -nj + wM 

= M a (l+ti + w^j + O (w 2 ) . (6) 

Combining eqs. (2) and (6) proves eq. (3). 

Hence the statistical strain U coincides with the elastic 
strain u. Note that in the (rare) cases where the micro- 
scopic displacements are affine, then £ = (I + w ) £q, \/£ 
This is a very strong statement since each link (and not 
only the average) obeys eq. (5); in this case, eq. (3) is 
easier to demonstrate. 

In summary, we have proposed a statistical characteri- 
zation of deformation. It averages the microscopic details 
of the current pattern, and of the reference state, to keep 
only the physical features statistically relevant at large 
scale. Hence different microscopic configurations which 
are statistically identical correspond to the same statis- 
tical strain. For instance, ductile metals like soft steel 
or aluminium have almost unchanged Young's modulus 
even much beyond their yield strain [19]; then, although 
different applied strains correspond to different micro- 
scopic structures, they have the same statistical strain, 
the same static stress and the same mechanical response 
for any physically reasonable static constitutive relation. 

Or, in a steadily flowing material, U is constant; we thus 
provide an operational definition for the thermodynamic 
stored strain that Porte et al. [20] introduced to allow 
a theoretical description of shear-induced phase transi- 
tions. 

This definition invites re-analysis of existing data, as 
well as experimental, numerical and theoretical tests. We 
have sucessfully performed such tests on an experiment 
which forces a two-dimensional foam to flow through a 
small constriction in a companion paper [21]. 
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